THE SPATIALLY HOMOGENEOUS RELATIVISTIC BOLTZMANN 
EQUATION WITH A HARD POTENTIAL 



HO LEE AND ALAN D. RENDALL 

Abstract. In this paper, we study spatially homogeneous solutions of the 
Boltzmann equation in special relativity and in Robertson- Walker spacetimes. 
£SJ i We obtain an analogue of the Povzner inequality in the relativistic case and use 

it to prove global existence theorems. We show that global solutions exist for a 
certain class of collision cross sections of the hard potential type in Minkowski 
space and in spatially flat Robertson- Walker spacetimes. 
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O ' 1. Introduction 

Under the assumption of homogeneity the non-relativistic Boltzmann equation 
has been extensively studied for many years. Homogeneity means that the unknown 
in the equation, which in general depends on time, spatial variables and velocity 
variables, is restricted to depend only on time and velocity variables. Since Carle- 
\ man gave the first proof of global existence [5] in the 1930s, many mathematicians 

^sD ' have obtained results on the homogeneous non-relativistic Boltzmann equation, and 

by now a rather complete mathematical theory is available. In the present paper we 
are interested in generalizations of this to special and general relativity. Compared 
to the non-relativistic case, the relativistic equations have not been studied much. 
Noutchegucmc and his colleagues obtained results on the homogeneous Boltzmann 
equation in [THJ [T9l [20] [21] in several different relativistic situations, but the scat- 
tering kernels they used are not physically well-motivated. The purpose of this 
paper is to obtain analogues of their results in more physically relevant cases. 

The scattering kernel is a quantity that determines the nature of collisions be- 
tween particles, and in the non-relativistic case several different types of scattering 
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kernel have been found to be of interest. For instance, the inverse power law gives 
the best-known types of scattering kernel, and they are further classified into hard 
and soft potential cases. In the relativistic setting it is not so clear which types of 
scattering kernel should be of interest, but a classification of (special) relativistic 
hard and soft potentials has been proposed in [SJ [25] by applying arguments sim- 
ilar to those used in the non-relativistic case. As in the non-relativistic case, the 
scattering kernels depend only on the relative momentum and scattering angle of 
two colliding particles. Consider a collision of two particles, and let p a and q a be 
their momenta before the collision, and p' a and q' a the momenta after the collision. 
Then the scattering kernel is given by a function of p a — q a and p' a — q' a . To be 
precise, the scattering kernel depends only on the following two quantities: 

(Pa ~ q a )( P a - q a ) and (p a - q a ){p' a - q' a ), 

where the indices are lowered by the Minkowski metric. The above two quantities 
are related to the relative momentum and the scattering angle respectively, and the 
relativistic hard and soft potentials are defined in terms of the above quantities. 

l 
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Their precise definitions will be given in Sections 11.11 and 11.21 On the other hand, 
the general relativistic Boltzmann equation does not seem to have been studied 
enough concerning specific types of scattering kernel. For instance, in [3J 2] the au- 
thors introduced a quantity S which is a function of x, p, q, p', and q' . The quantity 
S was used to play the role of the kernel of the collision operator of the Boltzmann 
equation, and strong assumptions on S were made to obtain local well-posedness 
of the Einstein-Boltzmann system. However, instead of introducing the abstract 
quantity, we may use the principle of general covariance to deal with the collision 
operator. If we write down the collision operator with an arbitrary Lorentzian met- 
ric replacing the Minkowski metric, then the special relativistic collision operator 
is extended to general relativity in a natural way. This argument is consistent with 
the following fact. There exists an alternative way to write the general relativistic 
Boltzmann equation. One can introduce an orthonormal frame on spacetime 
and parametrize the mass shell as p a = e"u^. Then, the collision operator no 
longer contains any explicit dependence on the metric and reduces to the special 
relativistic collision operator. We refer to [16] for a more detailed discussion of this 
approach. 

In this paper, we follow the procedure just mentioned. The general relativistic 
collision operator has the same form as in the special relativistic case, but a certain 
type of Lorentzian metric will replace the Minkowski metric. The definitions of hard 
and soft potentials will be understood by using the principle of general covariance, 
and eventually we will show that the Boltzmann equation has a global solution 
in the case of a certain type of hard potential in spatially flat Robertson- Walker 
spacetimes. We refer to [25l [31] for basic information about general relativity and 
to [3 [SI [TOJ and Chapter X of [7] for relativistic kinetic theory. For classical 
theories of the non-relativistic Boltzmann equation we refer to J6[ [12] and their 
references. 

1.1. Remarks on the relativistic Boltzmann equation. Since the Boltzmann 
equation is an equation describing the dynamics of collisional matter, it is necessary 
to understand the collision processes between particles in order to investigate the 
equation. The main difference between the non-relativistic and relativistic collision 
processes is that energy and momentum conservation in the non-relativistic case 
are replaced by an energy-momentum conservation, and this causes a difficulty 
in parametrizing post-collisional momenta in the relativistic case. Let p a and q a 
be two four-vectors describing momenta of two colliding particles in a relativistic 
situation, and suppose that the two particles produce momenta p' a and q' a after 
the collision. Then, energy-momentum conservation is written as 

p a +q a =p' a + q' a . 

Moreover, if all the particles are assumed to have the same mass, then the mass 
shell conditions 

are additionally imposed. Since the Greek index a runs from to 3, we have six 
constraints and two free parameters for the post-collisional momenta, and for the 
free parameters we use ui € § 2 as usual. Consequently, the post-collisional momenta 
p' a and q' a can be parametrized in terms of the pre-collisional momenta p a and q a 
with an additional parameter w G § 2 . There are several different ways known to 
parametrize post-collisional momenta in the special relativistic case, see for instance 
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[13l [29] . However, the parametrization suggested in [16] will be used in this paper 
because it can be applied effectively to the general relativistic case. Suppose that 
p a and q a are given, and dehne 

n a := p a + q a and t a := — n uj) 

for u€§ 2 , Note that t a is a general form of vectors orthogonal to n a , i.e. n a t a = 
for any to £ § 2 . The post-collisional momenta are represented as 

(1.1) p = ; and q = ■ 

[ } P 2 2 2^F 

It can easily be checked that they satisfy the mass shell condition and energy- 
momentum conservation. The scalar quantity g is called the relative momentum 
and is defined by (|1.3[) below. 

In this paper we are mainly interested in the Robertson- Walker spacetimes but 
first we consider the related case of Minkowski space. The spatially homogeneous 
Boltzmann equation in Minkowski space is written as follows: 

(1.2) d t f = Q(fJ) := f [ v$a(g,8)(fft - ff.)y/W\*"dq, := ^j, 

where /' = f{t,p'), ft = f(t,q'), f = f(t,p), and ft = f{t,q). Here, / is called 
the distribution function, Q the collision operator, the M0ller velocity, a the 
scattering kernel, 6 the scattering angle, and rj the determinant of the Minkowski 
metric. The first component p° can be solved for in terms of the other components 
p l due to the mass shell condition as 



where p = (p 1 , p 2 , p 3 ) denotes the spatial projection of p a , and the scalar quantities 
s and g are given by 



(1.3) s:=-( Pa + q a ){p a + q a ), g := ^ {p a - q a ){p a - q a ). 

They are called the total energy and the relative momentum respectively. Note 
that they are conserved quantities in the collision process, i.e. 



s = -(Pa + q' a )(p' a + q' a ), g = V(p'a'q' a )(p' a -q' a )- 

The scattering angle 9 is finally defined as the angle between p a — q a and p' a — q' a . 
It can be expressed in terms of Lorentz invariant quantities as follows: 

(1.4) cos* := = 1 _ afr.-iW-P"), 

9 s-4 

where s and (p a — p' a ){p a — p' a ) are two of the Mandelstam variables. Note that 
the scattering angle does not depend on how the post-collisional momenta are 
parametrized. 

In the second part of the paper we will consider the case of the spatially flat 
Robertson- Walker spacetimes. For a given scale factor R(t), the Boltzmann equa- 
tion is written as 



dtf~2-}^p l d p ,f^Q(fJ), 

i—l 



where R denotes the time derivative of R. The collision operator Q is the same as 
in (|1.2|) . but all the quantities in it such as (|1.3|) and (|1.4|) are calculated through 
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the Robertson- Walker metric. For instance, \f\q\ is replaced by R 3 (t), and indices 
are lowered by the Robertson- Walker metric, so we have 

p° = v/i + i? 2 (*)b| 2 

by the mass shell condition. However, in this paper we will take a different ap- 
proach. It is well-known that in the Robertson- Walker case the Vlasov equation 
can be explicitly solved for any given initial data. If the distribution function is 
expressed in terms of covariant variables, then it is independent of time, hence we 
have f{t,pi) = fo(Pi) hi the Vlasov case. Similarly, we can use covariant variables 
for the Boltzmann equation. Let us consider again the collision process by using 
covariant variables. To make the difference between contravariant and covariant 
variables clear, we will use v a for covariant variables, i.e. 

v a :=p a and v = v 2 , v 3 ). 

Then, the energy-momentum conservation is written as 

v a + u a =v' a + u' a , 

and the post-collisional momenta v' a and u' a are given by 

/ v a -(- u a g t a . v a -)- u a g t a 

zi = , and u„ — ., . 

2 2^W 2 2^F' 

where g and t a are basically the same as in the contravariant case, but they should 
be understood as quantities which are constructed from v a and u a . Consequently, 
if the distribution function is regarded as a function of t and v, instead of p, then 
the Boltzmann equation is written as 

(1.5) d t f = Q(f,f) = R- 3 [ [ V4<r(g,9)(ffl- ff*)dudu, v+ = 

where /' = f(t,v'), fl = f(t,u'), f = f{t,v), and /* = f{t,u). The first component 
v o is solved for in terms of the other components using the mass shell condition as 

v = -v/i + i?- 2 (t)M 2 , 

and the scalar quantities s and g and the scattering angle 9 are the same as (|1.3[) 
and (|1.4p respectively. The factor R~ 3 comes from the relation 

R 3 (t) dq = R- 3 (t)du. 

In the present paper the Minkowski case and the Robertson- Walker case will be 
studied separately. In the former case the Boltzmann equation will refer to (|1.2[) . 
while in the latter case the equation will refer to (|1.5|l . Comparing the two forms 
(jl.2[) and (|1.5[) , we might expect that results in the Minkowski case could be easily 
extended to the Robertson- Walker case, and this will be done in Section [3] 

1.2. Assumptions of the paper. The scattering kernel a is a function of the 
relative momentum g and the scattering angle 9. According to the way that the 
scattering kernel depends on its variables, it is classified into hard and soft po- 
tentials. The classification of hard and soft potentials in the relativistic case was 
originally introduced by Dudynski and Ekiel- Jezewska [3] and recently reformulated 
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by Strain [28] as follows: for soft potentials we assume that there exist 7 > — 2 and 
< b < min{4, 7 + 4} satisfying 

—p ) 9 , ~ 6 °'o(w) < a(g,u>) < g^ b a (uj) 7 



is 



<Jq(lo) < sin 7 1 



while for hard potentials we assume that there exist 7>— 2,0<a<7 + 2, and 
< b < min{4, 7 + 4} satisfying 

g a a (u)) < a(g,u}) < {g a + g- h ) a (cj), 
<To(ui) < sin 7 6. 

Here for any two quantities A, B the relation A < B mens that there exists a 
constant C such that A < CB. 

The scattering kernel. In this paper we assume that the scattering kernel has 
the form 

(1.6) cr(g,u) = 5 a sin 7 (9, -2 < 7 < -1, < a < 7 + 2. 

Since (^y) is a bounded quantity (note that s = 4 + g 2 ), a scattering kernel of this 
form falls into the hard potential case. 

Throughout the paper we use weighted L 1 spaces in momentum variables. Let 
i 1 (M 3 ) be the usual Lebesgue space of integrable functions on M 3 . Ll(R 3 ) denotes 
the weighted L 1 space with norm 



(1-7) \\fh,r~ / f(y)(y) r dy, (y) ~ VT+W 



The notations which will be used in this paper are as follows: the Minkowski metric 
is given by diag(— 1, 1, 1, 1), and the Robertson- Walker metric is 

ds 2 = -dt 2 + ^{^((dx 1 ) 2 + {dx 2 ) 2 + (dx 3 ) 2 ), 

where the scale factor R is a function of time t satisfying 

R{t) > 0. 

For simplicity we assume 

R{0) = 1. 

Greek indices run from to 3, and Latin indices from 1 to 3. p a denotes a four 
dimensional vector, while p denotes a three dimensional (contravariant) vector. 
Similarly, v a denotes a four dimensional vector, while v denotes a three dimensional 
covariant vector. To be precise, 

P = {p 1 ,? 2 ,? 3 ), v = (vi,v 2 ,v 3 ), 

while 



p° = ^l + R 2 (t)\p\ 2 , v° = -v = Vl + R- 2 (.t)\v\' 

Note that 



p° - «°. 
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The Einstein summation convention is used as p a p a = ^2 a =oP a P a , where the 
indices are lowered by p a = g a pp^ for a metric g a p . In some places we use pip 1 — 
Sj=i PiP 1 ■ The usual inner product • will only be used for three-dimensional vectors: 

3 3 

p-q = ^p l q\ v ■u = ^2v l u l , 

i=l i=l 
and | • | is such that \p\ 2 — p ■ p and \v\ 2 — v ■ v as usual. 

2. Existence results in the cases of bounded kernels 

In this part we briefly review the results of [181 [T9l [20l [21] , where the kernels of 
the collision operators are assumed to be bounded. For instance, the homogeneous 
Boltzmann equation in Minkowski space is written in |21j as 

(2.1) dtf = p^JI S{p ' q > p '> q 'Wf* - ?M du % 

An unknown quantity S is introduced to play a role of a kernel of the collision 
operator. The collision kernel S is assumed to be bounded uniformly on pre- and 
post-collisional momenta with an additional symmetry assumption: 

(2.2) < S < C, 

(2.3) S(p,q,p' \q') = S(p',q',p,q), 

where C is a fixed constant. The solution obtained in [3T] should be understood as 
a mild solution. By integrating the equation (12.11) with respect to the time variable 
from to t, we obtain the following integral equation: 

(2.4) f(t,p) = fo( P )+ f \ II S( P> q,p',q')(ffi-fU)dw^ds. 

Jo P J J 1 

By saying that / is a solution to the Boltzmann equation (I2.1j) , we mean / satisfies 
the integral equation (|2.4[) in L 1 -sense for each t. The main theorem can be stated 
as follows. 

Theorem 2.1 (Noutchegueme and Tetsadjio, [2T ). Suppose that the collision ker- 
nel S satisfies (|2.2p - (|2.3l) for a given positive constant C. Let r 6 (0, 56 ^. c ] and 
fo G X r be given for 

X r := L\R 3 ) n {/ > 0, a.e.} H {\\f\\ L i < r}. 

Then, the Cauchy problem for the homogeneous Boltzmann equation (|2.4p in Minkowski 
space has a unique global solution f G C([Q, oo);X r ) satisfying 

sup ||/(t)|Ui < H/olUi. 

t£[0,oo) 

The above result is extended to the cases of curved spacetimes in [19] 20j , and the 
initial condition is improved to any arbitrary large initial data (note that the initial 
condition in the above theorem depends on the constant C, which is the upper 
bound of the collision kernel). The equation considered in 20 can be written as 

(2.5) d t /-2|5> i; cV/=^ J! A(t,p,q,p',q')(f'f:- ffjfa^, 
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where the scale factor R = R(t) is assumed to be given such that it is differentiable 
and bounded from below, and the collision kernel A is similarly assumed to be 
bounded with a symmetry assumption: 

(2.6) R(t) > C, 

(2.7) < A < C, 

(2.8) A(t,p,q,p',q')=A(t,p',q',p,q). 



The solution of the equation ([2. 5ft obtained in [20] should also be understood as a 
mild solution. We define a characteristic curve of (|2.5|) by 



(2.9) ^ = _2 | p with p(o) = v. hence p(*) = R~ 2 (t)y- 



Then, the equation (|2.5[) can be written as 

— /(t,p(i)) = (collision term), 
dt 

and by integrating the above equation in a similar way, we obtain 
f(t,p(t)) = fo(y) + / (collision term) ds. 



However, if we regard / as a function of t and y, where y and p = p(t) are related 
to each other by (|2.9[) , i.e. 

~r/(*)J/) — (collision term), 
then we obtain a different mild form, 



/(*, y) = /o(y) + / (collision term) ds, 
Jo 



and this was the argument of [5D]. To be precise, we express the above integral 
equation in terms of y and z ■= R 2 (s)q, and obtain the following mild form of the 
Boltzmann equation in the Robertson- Walker spacetime: 

(2.10) f(t,y) = fo(v)+J J J A(s,y,z,y',z')(ffi-ff,)dLJ^^ds. 

The solutions obtained in [20] are functions which satisfy the integral equation 
(|2.10|) in I/ 1 -sense for each t. If we recall that in this paper the Boltzmann equation 
in the Robertson- Walker spacetime refers to (|1.5j) and that to derive this form 
of the equation we considered covariant variables v instead of p, then we can see 
that considering the characteristic curve corresponds to considering the covariant 
variables and that the new variable y of (I2.10P corresponds to the covariant variable 
v of (ll.5p . The main theorem of [20] can be stated as follows. 

Theorem 2.2 (Noutchegueme and Takou, [IS]). Suppose that the collision kernel 
A satisfies (|2.6p ~ (|2.8[) , and let f £ L^IR 3 ) with fo>0be given. Then, the Cauchy 
problem for the relativistic Boltzmann equation in the Robertson-Walker spacetime 
(12.101) has a unique global solution f £ C([0, oo); L 1 (]R 3 )) with f(t) > satisfying 

sup ||/(t)|| L l < H/olUi. 

te[o,oo) 
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In the above result the scale factor R(t) is only assumed to be bounded from 
below, hence Theorem l2.2l includes Theorem l2.ll A similar result has been obtained 
in [19] in a little more general spacetime, which is a Bianchi Type I spacetime, and 
finally the Einstein-Boltzmann system has been studied in |18J . The solution space 
used in [18] is a weighted L 1 space, which corresponds to L}(R 3 ) in the notation 
of the present paper. If we combine the arguments and proofs of [HI HH [20], the 
following theorem is obtained. 

Proposition 2.1. Consider the Boltzmann equation (|2.10[) in the Robertson-Walker 
spacetime. Suppose that the collision kernel A satisfies (|2.6p - (|2.8p . and let /o G 
i}(M 3 ) with fa > be given. Then, the Boltzmann equation has a unique global 
solution f e C([0, oo); i}(R 3 )) with f(t) > satisfying 

Sup ||/(t)|| L ; < 11/olLj. 

t£[0,oo) 

3. Existence results for the hard potential cases 

In this section the existence result of Proposition 12.11 will be extended to the 
hard potential case. Proposition 12 . II shows that the Boltzmann equation (|2.10l) has 
a solution in the weighted function space L}(M 3 ) with the weight function (y). By 
assuming that the collision kernel A is independent of t and the scale factor satisfies 
R(t) = 1, we obtain the existence result in L}(K 3 ) for the equation (|2.4j) . This will 
be extended to the hard potential case in Section 13.21 The existence result for the 
Robertson- Walker case will be studied in Section 13.31 

3.1. Preliminaries. In this section we collect several lemmas. 

Lemma 3.1. The following inequalities hold in the Robertson-Walker spacetime: 



(3.1) 9 < yfs, 

(3.2) g < 2VpV, 

(3.3) < 2vW, 

(3.4) g<R(t)\p-q\. 



Proof. Since s = g 2 + 4, we obtain the first inequality. The second inequality is 
given by the first and the third ones. For the third inequality we notice that 

s = 2 + 2p°q a -2R 2 (t)(p- q) 

= 2p Q q° + 2y/l -2R 2 {t){p-q) + R i {t){p-qY 

< 2p Q q° + 2^1 + R2(t)(\p\ 2 + \q\ 2 ) + R*(t)\p\ 2 \q\ 2 

= Ap Q q\ 

and this proves the third inequality. The last inequality is also clear because 

9 2 = "(P° - q°) 2 + R 2 (t)\ P - q? < R 2 (t)\p - q\ 2 , 
and this completes the proof. □ 

Lemma 3.2. For any uj e § 2 , we have t a t a > R 2 (t)s in the Robertson-Walker 
spacetime. 
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Proof. In the Robertson- Walker spacetime we have 

t a = ( ni uj\ -n uj) = (R 2 {t){n ■ uj),n°uj), 
and the proof is a direct calculation: 

t a t a = -R 4 (t)(n ■ uj) 2 + R 2 (t)(n ) 2 > -R 4 (t)\n\ 2 + R 2 (t)(n ) 2 
= R 2 (t)((n ) 2 - R 2 {t)\n\ 2 ) = R 2 (t)s. 
This completes the proof. □ 

The following lemma is a well-known fact in the special relativistic case |13l , and 
we show that the lemma holds in general relativistic cases also. The lemma is proved 
by the same argument as in [13] . but we present it for the reader's convenience. 

Lemma 3.3. Let (j>' a ,q' a ) and (j> a ,q a ) be pre- and post-collisional momenta re- 
spectively, and consider the collision map (p a ,q a ) — > (p /a ,q /a ). Then, the Jacobian 
is given by 

d( P ',q') _ pWo_ 

d(p,q) Poqo' 

Proof. To prove the lemma, we use a different parametrization from (jl.ll) . 

p' a =p a +2^-t Q , q' a = q a - 2^-t a , 
where t a is the same as in Ijl.ljl . For convenience we write 

p' k =p k +Acj k , q lk = q k -Au\ A = -2^no. 

By the same calculations as in [T3] we obtain 

d(p', q') _ ( 8) + (d p3 A)u>* (d q] A)u>* 



d(p,q) " V ~(dpiA)^ 8)-{d qj A)w* 
= det (8) + {d p] A-d q] A)u l ) 
(3.5) = l + (d p ,A-d q *A)Lj\ 

We differentiate the conserved energy 



p m + q' Q =p° + q 



with respect to p 3 , and multiply lo j to obtain 

p' k dp' k q J k dq Jk \ )J _ P 4 



Po <V q'a dp 3 J Po 



LO = -U 



Po PoJ \Po %' 



Similarly we obtain 



Hence, (|3.5[) is given by 



( Si 


4V- 


(Pi 




\qo 




U 





d(p,q) \p' q' J \Po qa 



(3.6) ^4 = ( 4^ - 4- fe ) 1 ( ^ - ^ 
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Recall that n a = p a + q a and t a = (rijOji , —nou) for w£§ 2 , and then the above 
quantities are written as follows: 

(3.7) — uji — = — — (<7o% — nofljW = —^—q a t a . 

po go poqo ' poqo 

Similarly we obtain 

v -r^ k - ^ k = -^WoPk + iWk - vWk - pWk)" k 

Po % Pa% 

(3-8) = S-T{q' n k - n a q' k )u k = -^q'j" = —^q a t a , 

Po% Po% Po% 

where we used the energy-momentum conservation and the following: 
t a q — t a q 2— ~y^t a t — taQ ■ 

We plug (|3.7p and (|3.8[) into (|3.6I) . and this completes the proof. □ 

Lemma 3.4. For the collision operator the following property holds in the Robertson- 
Walker spacetime: for any measurable function k depending only on g, s, and U), 
we have 

^^(f'f:-ff*)(p Q Yd^dqdp 
p v q" 

\ ill ^^//•((^) P + (^) r -^) P -(^) P ) dw **- 



Proof. We use Lemma 13.31 to make the change of variables between pre- and post- 
collisional momenta as follows: 

(3.9) -^-rdpdq = -——dp'dq', 

and note that g and s are invariant quantities under the collision process and 
symmetric for p and q. Hence, the gain term can be written as 

l^.ff^ rdudqdp= HI k -^rap°Yd,dq>dp> 

k{g,s,uj) , 

n-n— //*(P ) dujdqdp. 

pUqu 

By interchanging p and q, it can also be rewritten as 

ti^lff Ap oy diodqdp = m h^l fWOydLjdqdp _ 

Hence, we obtain the following representation for the gain term: 

ti^lff Up oy dujdqdp = l _m M^Uff^y + iq >° r)dujdqdp . 

After applying the same argument to the loss term, we obtain the desired result. □ 

Lemma 3.5. Consider the collision process in the Robertson- Walker spacetime. Let 
(p' a ,q' a ) and (p a ,q a ) be pre- and post-collisional momenta respectively. Consider 
the following quantity for r > 1; 

G = (p'°Y + (q'°) r - (p°) r - ( q °y. 
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Then, G satisfies 

(3.10) g < a((p r-y +p°(q°y- 1 ). 

If oj is restricted to a subset {oj € S 2 : \n ■ cj\ < "^75 ML then G satisfies 

(3.11) G < C r ((p ) r -i(q )i + (p )Hq r-^ - c r ((p°Y + (q°) r ), 

where C r and c r are two different positive constants depending on r. 

Proof. Note that p'° + q' a = p° + q° is a conserved quantity for each p° and q° . Let 
p a and q a be given. Then, G reduces to a function of u> and can be written as 

G(w) = (p ,0 ) r + (q ,0 ) r - (p°Y - (q°) r 

= (p' a ) r + (p + q -p' Y-(p Y-(q ) r - 

Note that G has minimum at p'° = (p° + q°)/2, i.e., when p'° = q' Q , and is mono- 
tonically increasing as p'° tends to or p° + q°. Hence, G attains its maximum 
when p'° — q'° is extremal. Without loss of generality we may assume p'° > q'°. 
From (ll.ip we obtain 

,o ,o t° R 2 {n-u) 

p -q = 9^== = 9- 



VtJ^ y/R?(n°) 2 - R 4 {n-uiY 

This quantity attains its maximum when oj is parallel to n, which means that p is 
largest and q is smallest when n ■ u> — \n\, in particular q 1 < min{p°, q }, which 
implies again (q'°Y < min{(p ) r , (g°) r }. Consequently, G attains its maximum 
when n ■ ui = \n\ and is estimated as 

G(u) < G{n/\n\) < (p'°) r L =n/|n| - max{(/) r , (q°Y} 



p°+q° g R 2 \n\ 



- max{(p°r, (q°Y}- 

aj—n/\n\ 



2 2^/tJ^ 
By applying Lemma 13. 11 Lemma 13.21 and the inequality 
(3.12) (a + bY < a r + b r + C r (a r - 1 b + ab r - 1 ) for r > 1, 



we obtain 



GM<(^ + ^±il) -max { (/r,( ? T} 

< (p + g ) r -max{(p ) r ,(q )'-} 

< min{(p°) r , (q )''} + C r ((p°) r -y +p ( g ) r - 1 ) 



and this proves the first result. 
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To prove the second result, we take the assumption \n ■ u>\ < ^75 M an d suppose 
p'° > q' Q . Then, p'° is estimated as 

1 R 2 \n-uj\ 



p'° < P +<1 



2 2 ^i?2( n 0)2_ i? 4( n . w )2 

< p° + q | g R 2 \p + q\ ^ p° + q | g 

2 2 ^2 ^i?2(p0 + g 0)2 _ 1^4^ + 9 |2 " 2 2 

< (y^+y^) 2 

2 

where we used Lemma I5TT1 Then, G is estimated as 

g < 2( P '°y - (p°y - (q°y 

< { ^^ r - (p°y - ( q °y 

< + + C^rW)* + (p°)*(fl°r *) - (p°) r - (l°Y 

< C r ((p°y-i(q°)i + (p°)'(g ) r -') - c,((?/T + (?°) r ), 

where we used (|3.12[) . The constants C r and c r are two different positive constants 
depending on r, and this completes the proof. □ 

Remark 3.1. Those types of inequalities given in Lemma \3.5\ are called the Povzner 
inequality, which was originally proved by Povzner in |23j . This inequality has 
been crucially used to prove existence theorems for the non-relativistic spatially 
homogeneous Boltzmann equation by Elmroth [11] and Mischler and Wennberg [17) . 
The sharpest form of the Povzner inequality is given by Mischler and Wennberg, 
but Lemma \3.5\ corresponds to a relativistic extension of Elmroth 's result. 



3.2. Hard potential case in Minkowski space. Comparing the two forms (|1.2|) 
and (|2.ip for the Boltzmann equation, we can see that the collision kernel S of (|2.1j) 
corresponds to gyfso~(g,6) of (|1.2[) . Hence, if the quantity g^fso~{g,6) is suitably 
truncated, then the truncated equation has a global solution by Proposition 12.11 
For simplicity the following notations will be used: 

min^v^m} . 
w*,m := n~n ' 9m ■= mm j,ra}, cro,m(w) := mm{a {u}),m}. 

For each integer m, let f m be a solution of the following truncated equation with 
initial data / m (0) = /o: 

(3.13) dtfm = Qm(fm, fm), 

where Q m is defined as 



Q m (h,h) '=11 «<£,m(Sro) o 0Q,mM(/i h, - hh*) du; dq. 

Then, the truncated equation has a unique global solution f m € C([0,oo); L\( 
by Proposition 12.11 Consequently, a sequence {f m } is obtained, and existence of 
solutions for the original equation (jl.2p is proved by showing that the sequence {/ m } 
is a Cauchy sequence in L}(]R 3 ). This argument was used for the non-relativistic 
case by Mischler and Wennberg in [17] , and below we will show that their argument 
is applicable to the relativistic case. 
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Lemma 3.6. For any r > and T > 0, there exists a constant C r which does not 
depend on m such that if ||/o||i,r is bounded, then 

Sup sup ||/ m (t)||l, r < C r . 
m te[0,T] 

Proof. We first note that by Proposition 12. II 

sup ||/ TO (t)||i, P < C for 0<r<l, 
te[o,oo) 

where C does not depend on m, and for r < s, 

||/m(*)||l,r < ||/m(*)||l,a. 

We now assume r > 1. By direct calculations we have 

j7||/m(£)||l,r 

= J J J ^,m(9m) 8 ffo, m (w)(44, - fmfm*){p°Y dudqdp 

= \JJJ v<pM9m) a o-o,Mfmf m *((p' a ) r + (<z'°) r - (P°Y - (q°) r )du;dqdp, 
where we used Lemma 13.41 We apply Lemma 13.51 to obtain 

J t \\fm(t)\\l,r<h+h-h, 

where 

h = C r [[[ V^9rn) a <JoM^)fnJ m *((p ) r ^Q +P°(q ) r ^)diodqdp, 

JJJ\n- U \>-^\n\ 

h = C r ({{ V4,, m (9m) a <Jo,Mf m frn*((p a ) r -Hq°^ + (p^Hq )^) dm dq dp, 

JJJ\nM<^\n\ 

h = C r / / / ^,m(ffm) O CT0,m(w)/m/m*((p ) r + (<7°) r ) dudqdp. 

J J J\n-u\<^\n\ 

The second term 1% is easily estimated by using Lemma 13.11 as 

h<C r JJ fmf m *(p°) r -^{q )^ dqdp < C r \\ f m (t)\\ 1>r _ i +f \\f m (t) || 1>|+f . 

Consider now ero,m( w )> which is defined by 

<7o. m (u;) := minjsin 7 9,m} for — 2 < 7 < — 1. 

Note that cro im (w) is integrable on § 2 for 7 > —2, and there exists a constant C 7 
satisfying L 2 o"o im (u;) < C 7 , where the constant C 7 does not depend on m. On 
the other hand, since 7 is negative, we have o~o im (uj) > 1 for any m. Moreover, the 
integration domain of I3 is a set with Lebesgue measure 

H{\n ■ u\ < \n\/V2} = 2\/2tt, 

which does not depend on m. Hence, I± and ^3 can be estimated as 

h<D r II v 4> , m (g m ) a f m f m *(p ) r ~ 1 <i dqdp, 
J Jm. 6 

h>d r {9m) a fmfm*(q°) r dqdp, 
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for some constants D r and d r . We now fix the constants D r and d r to split the 
domain by {£) r (p°) r - 1 < rfr(<?°) r_1 } and {£) r (p°) r - 1 > rf r (9 ) r_1 }, and then obtain 
/i < /n + /12, where 

J n = D r v 4> , m (g m ) a f m f m , t (p y- 1 q dqdp, 

J J D r ( P °y-i<d r (q°y-i 

hi = D r \\ v 4>tm (g m ) a f m f m 4p°) r ^ 1 q°dqdp. 

J J D r (p°r- i >dAq°r- 1 

We now obtain 

hi < h- 

In the case of /12, we may simply use (j m ) a < C(p°q°)% < C r (p a ) a - Then, J12 is 
easily estimated as 

Il2 < Cr||/m(i)|jl,r-l+a- 

Combining the above estimates, we obtain 

|ll.U*)lkr 



< 



C, (||/ ro (i)lll, r -I + f H/-Wlll,i + f + ll/rr.Wlll.r-1+a) < C r \\f m (t)\\ hr , 



where we used the fact that < a < 7 + 2 and — 2 < 7 < — 1. Then, the lemma is 
proved by applying Gronwall's inequality. □ 

Lemma 3.7. Consider the sequence {f m } on any finite time interval [0,T]. For 
any small number 5 > 0, there exists a positive integer M such that if k,m > M , 
then 

sup ||/ fc (i)-/„,(i)||i,i< J. 
te[o,T] 

Proof. Let k < m be two positive integers. By direct calculations, 

^ll/fcW-/m(*)||l,l=/ + J, 

where 

Sgn(/fe - f m ){Qk(fk,fk) - Qk{fm,fm))P° dp, 
J = j Sgn(/fe - f m ){Qk(fm, fm) ~ Qm(fm, fm))P° dp, 

and / and J will be estimated separately. The first term / is split again as 

1 = \ f ff S S n tfk - /m)w0,fc(3fc) a Cro,fe(w) 



2 . . 

x {(fk — fm)(fk* + fm*) + {fk + fm)(fk* ~ fm*) 

- {fk - fm)(fk* + fm*) - {.fk + fm)(fk* ~ fm*))P° doJ dq dp 

=: h + I 2 + h + h. 
Each Ii is estimated as follows: 

h<lJJI v^M^oMlf'k - f'mKf'k* + fL)P°du)dqdp 
= 7; j jj ' v<p, k (gk) a o-o,k(u)\fk - fm\(fk* + fm*)p'° dujdqdp, 
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h^lfff ^k{9k) a voM{fk + O/L - fL*\p°dudqdp 
= ^ JJJ v^ k (g k ) a a , k (u})(f kif + f m *)\fk - fm\q'° dujdqdp, 

h = ~ j Jj sgn(/ fe - fm)v4,,k{gk) a vo,k{u)(f k - f m )(fk* + f m *)p° divdqdp 

= ~ JJJ v<p,k(gk) a vo,k(u)\fk - fm\(fk* + fm*)p° dcvdqdp, 
and finally 

h < ^ J J J «0,fc(3fc) a cro,fe(w)(/fe + f m )\fk* - fm*\P° dwdqdp 

< ^ J J J v<t>,k(9k) a (ro,k( w )(fk* + fm*)\fk - fm\q° dudqdp. 

Therefore, I is estimated as 
1 
2 



I< \ III Vt,k{gk) a G0,k{u)\fk - fm\(fk* + fm«)(p /0 + q'° ~ P° + q") dudqdp 



<C'JJ (g k ) a \fk - U(fk* + fm*)q° dqdp, 



where we used 
By using g k < 2y / p°<Z°, we obtain for / 



p' + q m =p + q °. 



I<C J J \fk-f m \(fk*+f m *)(p°)%(q°) 1+ %dqdp 
(3.14) <Crsup||/ n (t)|| lil+ .||/ fc (t)-/ m (t)|| 1 ,.. 

n 

To estimate the second term J, we note that 

|^,fe-^,m| = ^^|min{5Vs,fc}-min{gVs,m}| 
minjgy^, m} 

< 1 {g^S>k} ^5 = 1 {g^>k} v 4>,r 

and similarly 

\(g k ) a - (g m ) a \ < i {g >k } (9m) a , 

|cO,fc(w) — CTo, m (w)| < l{ s i n 7 0>fc}<7'O,m (<*>). 

Hence, J can be estimated as 



J < J J J |w0,fc(5fe) a CT O ,fe(w) - V0,m(5m) a oo,mM||/^/^„ - f m f m *\p°dwdqdp 

-ill 1 {gV^>k} V 4>,m(gk) a ^0,k( ljJ )(fU'm*+ fmfm*)p° du} dq dp 
*-{g>k}V4>, m {9m) a ^0,k{^){f' m f m * + fmfm*)P° doo dq dp 



+ 



III 1 ^ 6 ^ v ^ m ^ 9m ^ aa ^ m ^^' m ^' m * + fmfm*)P° dwdqdp 
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Note that each J; can be separated into two terms: a gain term containing f' m f' m . Jf 
and a loss term containing f m fm*- The gain and loss terms are estimated in the 
same way after making the change of variables (p, q) ■<-> (p',q'), hence we only 
present the estimates for the loss terms. To estimate Ji, we take a small number 
e > and use g\fs < 4p°q° from Lemma 13. II 



Jl<C J J l{4 P °q°>k}(gk) a f m f m *P°dqdp 

<C // l{4p°g°>fc}/m(p°) 1+f jm*(g°) f dq dp 



<^JJ Mip^>k}fm(p ) 1+ ^ +£ fm*iq°^ +£ dqd P 

(3-15) < ^■||/ m (t)|| 1 , 1+§+£ ||/ m (f)||i, i+£ . 

To estimate J2, we use g < \p — q\ to obtain 

h<C J J l { \ p - g \>k}fm(p ) 1+§ f m *(q°)% dqdp 

~ C ll 1 {bl>#}u{kl>|}/™^°) 1+f /™*^ ) f dqdp 
< C||/ ro (*)||i,i+§ J l {k[ > k } f m *(q )z dq 
(3-16) < j ||/ m (i)||?, 1+f • 

For J3 term, we use sin6> m 8 for < 8 < ^. Hence, the condition sin 7 8 > k is 
equivalent to 8 < Chi since 7 is negative. We first estimate J3 as 

J 3 <CJJJ l {9<ck ± } a , m (")fm(P°) 1+§ f m *(q ^ dcodqdp 
< C||/ m (t)||l,l+|||/ m (t)||l,| y l {e<cfc l } f7 ,m(w)dw. 



The integration on § 2 above is estimated as 



/i>Ch~t 
1 1 a , m (cv) dw < 2tt / sin 7+ 1 6>d6» < CTc^, 

where the constant depends on 7. Note that —1 < (7 + 2)/7 < 0, and the third 
term J3 is estimated as 

(3-17) J 3 < Cfc^||/ m (i)||l,l+f H/mWIkf ■ 

We combine (|3.14j) . (|3.15j) . ()3.16[) . and ()3.17|) . and apply Lemma f3~6l on any finite 
time interval [0, T] to obtain 

j t \\h{t) ~ < C(Ar £ + fc" 1 + + C||/ fc (i) - / m («)||i, s 

< C(k- £ + k'-^) + C\\f k (t) - / m (t)||i,i. 

Since /fc(0) = / m (0) and (7 + 2)/j is negative, we obtain the desired result by 
applying Gronwall's inequality. □ 

We now obtain the following theorem. 



THE BOLTZMANN EQUATION 



17 



Theorem 3.1. Suppose that the scattering kernel has the form (|1.6|) and initial 
data satisfies fa G L^M 3 ) for r > l + § with /o > 0. Then, the Boltzmann equation 
(|1.2[) m Minkowski space has a unique global solution f € C([0, oo); L}(R 3 )) with 
f(t) > 0. 

Proof. Lemma f3 . 71 shows that the sequence {/ m } is a Cauchy sequence in L}(R 3 ). 
Hence, there exists a solution / to the Boltzmann equation (|1.2[) . The initial 
condition /o G L*(R 3 ) with r > 1 + § comes from (|3.15[) . and nonnegativity of the 
solution is guaranteed by Proposition 12.11 Uniqueness of solutions can be easily 
proved by following the calculations given in the proof of Lemma 13.71 The initial 
condition of the theorem is rather strong in the sense that L\ solutions are obtained 
from L\ initial data for r > 1 + hence the proof of uniqueness is much easier 
than that of [17]. This completes the proof of the theorem. □ 

3.3. Hard potential case in the Robertson- Walker spacetime. In this part 
we extend the existence result of Theorem 13. II to the Robertson- Walker case. The 
argument is basically the same as in the Minkowski case. We first truncate a certain 
part of the collision kernel, then existence for the truncated equation is guaranteed 
by Proposition [2Tj We obtain a sequence of solutions to the truncated equations, 
and by showing that the sequence is a Cauchy sequence we obtain a solution. 

To deal with the Boltzmann equation (ll.5[) in the Robertson- Walker spacetime, 
we use the covariant variable v. In the Minkowski case we have p° = (p), so there 
is no difference between the following quantities: 

J m P )(p°rdp= J f(t,p)( P ) r dp. 

However, in the Robertson- Walker case we have 

v° = ^l + R- 2 {t)\v\ 2 and (v) = + \v | 2 . 
If we introduce a new quantity: 

|/(*)kr := j f(t,v)(v°ydv, 
then we have the following relation: 

(3-18) |/(t)|l,r<||/(t)||l,r< fl r (*)|/(*)|l,r, 

where we used 

(3.19) v° < (v) < R(t)v°. 

Note that 

l/(0)|l,r = |l/(0)||l,r, 

since we assume R(0) = 1. 

The main goal of this section is to extend Theorem 13. II to the Robertson- Walker 
case. We will show that the lemmas in the previous section can be applied to the 
Robertson- Walker case. To do this we first use v° as a weight function, i.e. we 
estimate |/(t)|i, r , because all the calculations in the Minkowski case are naturally 
extended to the Robertson- Walker case when using v° instead of (v) , for instance we 
can use the energy-momentum conservation. Then, by the relation (|3.18l) we obtain 
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estimates for ||/(i)||i,r- Similarly to the Minkowski case we modify the Boltzmarm 
equation (|1.5[) as 



dtfm = Qm(fmjm) ■= R 3 J J Vj >m (g m T^O.rn (w) {f m fm* ~ fmf m *)dujdu, 

where 

mm{gy/s,m} . 
v<t>,m = 5-75 , Jm=min{g,m}, a , m = mw.{<T {u>), m}. 



Then, the truncated equation has a unique global solution by Proposition [2TTJ The 
following lemmas show that Lemma [5751 and 15771 can be extended to the Robertson- 
Walker case. 

Lemma 3.8. For any r > and T > 0, there exists a constant C r which does not 
depend on m such that if ||/o||i,r is bounded, then 

sup Sup \f m (t)\l,r + \\fm(t)\\l,r < C r . 

m te[Q,T] 

Proof. We first estimate \f m (t)\i,r> an d then obtain the desired result by using the 
relation (|3.18[) . By Proposition 12. II and the relation (|3.18[) we have 

sup \f m {t)\i, r <C for 0<r<l, 
teto.oo) 

and now assume r > 1. In the Robertson- Walker case, v° depends on time and 
decreases as time evolves for each v. To be precise, 

v° = v/l + i?- 2 «M 2 and dt v" = -^ l -^-<0, 

since we assume R(t) > 1 and R(t) > 0. If we follow the calculation of the proof of 
Lemma l3.6i then we obtain 



-Tj_\fm(t)\l,r 

= jjj V <t»*n{g m ) a CTO, m (u)f m f m *((v' y + (U ,Q Y - (V°Y - (U a ) r ) dujdudv 

+ J f m {t,v)^- t [{v°y] dv, 
and the second integral is negative. Hence, we may only consider 

jll/m(*)|l,r 

< JJJ v^ m (g m ) a a , m (u)f m f m 4(v'°y + (u ,Q y - (v°y - (u y)dwdudv, 

and follow the same calculations of Lemma 13.61 to obtain 

sup sup |/ m (t)|i, r < C r . 

m te[0,T] 

Consequently, the relation p.!8[) gives the desired result. □ 
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Lemma 3.9. Consider the sequence {f m } on any finite time interval [0, T], For 
any small number S > 0, there exists a positive number M such that if k,m > M , 
then 

sup ||/jfc(i)-/ m (t)||i,i<*. 
te[o,T] 

Proof. Similarly to the previous lemma, we first estimate \fk — /m|i,i and then 
obtain the desired result by using (|3.18[) . By direct calculation, we have 

^|/fcW-/mWU,l = ^ / \fk(t,v) - f m (t,v)\v° dv 

= fd t \\f k (t,v) - f m (t,v)\\v° + \f k {t,v) - f m (t,v)\d t v°dv 



Sgn(/fe - fm){Qk{fk, fk) ~ Qm(fm, fm))v° dv 



&(t) 



\fk(t,v)-f m (t,v)\ l -+dv 



< / Sgn(/ fc - fm){Qk(fk, fk) - Qm(fm, fm))v° dv. 



Hence, we can follow the proof of Lemma 13.71 and obtain a positive number M such 
that if k,m > M, then 

SUp \f k (t) - /m(*)|l,l < S. 

te[o,T] 

Consequently, we obtain the desired result by applying p.!8p . and this completes 
the proof. □ 

By the same argument as in Theorem 13. 11 we obtain the following theorem. 

Theorem 3.2. Suppose that the scattering kernel has the form of (|1.6[) and initial 
data satisfies fo £ Lj;(R 3 ) for r > 1 + § with fo > 0. Then, the Boltzmann 
equation (|1.5[) in the Robertson- Walker spacetime has a unique global solution f £ 
C([0,oo);L\(R 3 )) with f(t) > 0. 

4. Summary and outlook 

In this paper global existence theorems have been proved for spatially homo- 
geneous solutions of the Boltzmann equation in Minkowski space and in spatially 
flat Robertson- Walker spacetimes. This was done for a class of collision kernels of 
hard potential type and these theorems extend existing results for smooth collision 
kernels with compact support to cases which are closer to those which naturally 
arise in physical problems. 

There are many directions in which this work might be generalized. Do the 
special relativistic solutions constructed here converge to equilibrium as t — > oo? 
If so, can a useful analogue be proved for the Robertson- Walker spacetimes? Note 
that in the latter case equilibrium is impossible since the existence of an equilibrium 
solution of the Boltzmann equation in a spacetime implies the existence of a timelike 
Killing vector field (PUj, p. 1167) and vector fields of this kind do not exist in most 
Robertson- Walker spacetimes. It is nevertheless the case that in practise ideas 
related to equilibrium matter distributions are used in cosmology and so it should 
be possible to formulate some mathematical analogue of equilibrium solutions in 
an expanding cosmological model. 
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Now that a global existence result has been obtained in a class of homogeneous 
and isotropic cosmological models it is natural to ask whether a similar result can 
be proved in general homogeneous cosmological models which expand for ever. In 
this paper the device of writing the equation in terms of the covariant components 
of the momentum is used to simplify the equation. In fact this trivializes the kinetic 
part of the equation. The same trick would work in more general models of Bianchi 
type I but not in general homogeneous models. It is nevertheless the case that 
this transformation does produce some simplification in general and this has been 
exploited in the study of the late-time behaviour of the Einstein- Vlasov system in 
[22]. 

Another generalization is to look at the global existence question for homoge- 
neous solutions of the Einstein-Boltzmann system. A template for this could be 
provided by the results of this type for the Einstein equations coupled to other 
matter models proved in [23]. For the class of collision kernels considered here 
there are no local existence theorems available for inhomogeneous solutions of the 
Einstein-Boltzmann system. There are global existence theorems available for the 
Boltzmann equation with this type of kernel in special relativity with small or close 
to homogeneous initial data. The small data problem for the Einstein-Boltzmann 
system is out of reach at present since even the corresponding problem for the 
Einstein- Vlasov system has not been solved. On the other hand in the presence 
of a positive cosmological constant there are global existence results for solutions 
of the Einstein- Vlasov system evolving from data which are close to homogeneous 
[26]. 

It is well-known that the Boltzmann equation is not well-posed in the past time 
direction. In cosmology this equation is of interest for the very early universe and 
so it is natural to enquire if solutions of the Einstein-Boltzmann system can be 
constructed which extend all the way back to the big bang. It is not reasonable 
to construct these by evolving backwards in time and a possible alternative is to 
pose data at the singularity. A formal study of this type of procedure has been 
carried out in |30j but corresponding existence proofs have not yet been developed. 
In conclusion, the study of the Boltzmann equation in curved spacetimes and its 
coupling to the Einstein equations gives rise to a variety of challenging mathematical 
problems. 
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